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Improved  Element-Free  Galerkin  Method  for 
Electromagnetic  NDE  Model 
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Abstract.  This  paper  presents  an  improvement  of  the  new  Element-Free  Galerkin  method 
(EFG)  in  NDE  application.  By  using  the  higher  order  orthogonal  basis  function  instead  of  the 
conventional  polynomial  basis,  the  ill-condition  and  inversion  of  the  shape  function  are 
eliminated.  Therefore,  the  accuracy  of  the  results  is  improved.  The  preliminary 
one-dimensional  (1-D)  and  two-dimensional  (2-D)  examples  are  presented  to  demonstrate  the 
accuracy  and  efficiency  of  the  improvement. 

Keywords.  NDE,  meshless  method,  EFG  method,  orthogonal  basis 


Introduction 

Nondestructive  Evaluation  (NDE)  methods  are  widely  used  in  industry  to  control  and 
maintain  product  qualities.  The  availability  of  a  numerical  model  is  extremely  important  to 
visualize  the  field-flaw  interaction  in  NDE.  Recently,  meshless  methods,  known  as  a  new 
numerical  approach,  have  been  developed  to  eliminate  the  disadvantage  of  the  conventional 
Finite  Element  Method  (FEM).  Without  the  reliance  of  a  mesh,  meshless  methods  can  be 
implemented  to  model  complex  geometries,  such  as  tight  crack  and  dynamic  corrosions. 

Element-Free  Galerkin  method  (EFG)  is  one  of  meshless  methods  by  utilizing  moving  least 
square  (MLS)  approximation  and  Galerkin  formulation.  The  MLS  relies  on  the  following  three 
components:  a  weight  function,  a  basis  function,  and  a  set  of  position-dependent  coefficients. 

We  discovered  that  the  shape  function  is  ill-conditioned  when  using  high  order  polynomial 
basis.  This  will  prevent  the  implementation  of  high  order  polynomial  basis  for  higher  accuracy. 
To  solve  this  problem,  an  improved  formulation  of  the  EFG  method  using  orthogonal  basis 
function  is  proposed  in  this  paper.  The  condition  of  the  shape  function  is  improved  compared 
with  the  polynomial  basis  and  the  convergence  of  the  iteration  solver  is  faster  and  the  accuracy 
of  the  solution  is  increased. 

1.  Formulations 

In  the  EFG  method,  a  set  of  nodes  are  used  to  construct  the  discrete  system  of  equations 
approximating  the  solution.  And  a  set  of  quadrature  points  are  defined  to  compute  the  integrals 
over  the  solution  domain  in  Galerkin  procedure. 

The  MLS  approximation  can  be  written  as 
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(1) 


uh(x)  =  ^i®J(x)Uj 

j= 1 

where  n  is  the  total  number  of  nodes  and  the  shape  functions  <t>y  are  given  by 


=  Y,P'(x){B  '{x)C(x^‘j  =  pTb  'cj  •  (2) 

i= 0 


where  B(x)  =  PTW(x)P  ,  C(x)  =  PTW(x) ,  W  is  the  weight  function,  P  is  the  basis  function  and 
m  is  the  order  of  the  basis  function. 

The  choice  of  basis  function  is  essential  for  the  stability  of  the  numerical  results  because 
there  is  an  inversion  of  matrix  B.  The  ordinary  polynomial  basis  function  is  straight  for 
implementation  but  results  in  an  ill-condition  for  high-order,  which  makes  the  solution  not 
convergent. 

However,  an  orthogonal  basis  can  be  used  to  solve  this  problem.  The  definition  and 
formulations  of  the  orthogonal  basis  are  shown  below: 


Vector  p(x)  =[p1(x),...,p  (x),...,pm(x)]r  is  an  orthogonal  basis,  if  we  define  the  inner 


product  as 


(pi  (x),  pk (x)}  =  ^w'(x-  xt ) p , (x,  )pk  (x, )  = 


i= 1 


0  (l*k) 

C,±  0  (/  *  k) 


Therefore,  its  shape  matrix  becomes 

{Pi(x),Pi(x)) 

B(x)  = 

( Pm(x),Pm(x )) 

The  recursive  relationship  of  orthogonal  basis  is  given  by, 


(3) 


(4) 


Pi  (x)  =  1 

Pi(x)  =  {x-a2)p\  (x) 

Pk+ i  (x)  =  (x  - ak+]  )Pk  (x)  -  pkPk_l  (x) 


(5) 


where  k  =  1,2,...,  m.  The  parameters  a,  p  are  obtained  as, 

_  (y’HV), /?*(*)) 

k+1  {Pk(x),Pk(x)) 


where  k  =  0,1,...,  m- 1 .  Substituting  (4)-(6)  into  (2),  the  shape  functions  are  then  expressed  as, 
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(7) 


“  \v(  X  —  X  i)  P  :  (x 

<b,(x)  =  Y- -  7  Pj(x) 

M  (PjW'PjW) 

The  partial  differential  form  of  the  shape  function  is 


c/0/  (x) 


dx 


■E 

j= i 


dwix-Xj )  dp  j(x) 

m  - - - Pjix^pAx)  m  \v{x-Xl)pj(Xl) - - 

Zdx _ _ _  dx 


7=1 


Pj(x\Pj(x))  (pjix^Pjix)} 

dlpj(x\Pj(x)) 
wix-x^Pjix^Pjix)^ 


dx 


j(x),Pj(xjf 


(8) 


And  the  partial  differential  of  inner  product  form  is 


d(pj(x),pj(x )) 


^  r  l 

Xvv(x-x')^/(x')r 


7=1 


TV 

•Y 

7=1 


Jw(a:  -  jcz-  ) 
dx 


PjiXi)]2 


(9) 


These  equations  (3)-(9)  indicate  that  the  shape  functions  in  EFG  method  avoid  the 
inversion  of  shape  function  matrix  when  orthogonal  basis  is  used.  This  method  makes  the 
implementation  of  EFG  method  easier,  also  increases  the  accuracy  of  the  solution  by 
improving  order  of  the  orthogonal  basis. 


2.  Implementation 

2. 1  One-Dimensional  Electrostatic  Problem 

A  1-D  electrostatic  Poisson  problem  for  the  potential  function  u  is  described  as  follows: 

\v2u=-x  (Q:0<X<1) 

j*(0)  =  0  «’( 1)  =  0 


EFG  method  is  implemented  with  polynomial  basis  and  orthogonal  basis  of  various  orders 

for  comparison.  The  exact  solution  is  u(x)  =  —  x -  —  x3  for  comparison  with  the  numerical 

2  6 


solutions. 

Fig.  1  shows  the  exact  solution  and  the  numerical  solution  using  EFG  method  with  1st  and 
2nd  order  polynomial  basis.  We  discovered  that,  although  the  solution  with  2nd  order 
polynomial  basis  convergences,  it  starts  to  show  unstable. 
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Figure  1.  The  comparison  of  exact  and  EFG  solution  with  1st,  2nd  order  polynomial  basis. 

Fig.  2  shows  the  exact  solution  as  well  as  the  solution  using  EFG  method  with  3rd  order 
polynomial  basis.  The  numerical  solution  is  no  longer  convergent  because  of  the 
ill-conditioned  shape  function. 


rr\ 

Figure  2.  The  comparison  of  exact  and  EFG  solution  with  3  order  polynomial  basis. 


In  Fig.  3,  the  numerical  solution  using  EFG  method  with  orthogonal  basis  function  is 
presented.  The  order  of  the  basis  function  varies  from  1st  to  4th;  the  numerical  solutions  are 
stable  and  keep  convergent,  which  validate  the  formulation  of  the  orthogonal  basis  function. 


Figure  3.  The  comparison  of  exact  solution  with  various  order  orthogonal  basis. 
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2.2  Two-Dimensional  Eddy  Current  Problem 


The  2-D  eddy  current  modeling  geometry  is  shown  in  Fig.  4.  Assume  an  infinite-size  2-D 
conducting  plate  with  a  tiny  crack  along  y  direction  (length  =  1)  is  placed  under  a 
time -varying  harmonic  current  source  in  the  x  direction. 
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Figure  4.  Modeling  geometry  for  2-D  problem. 

The  induced  current  will  be  along  -x  direction  in  the  absence  of  crack,  and  will  be 
bended  by  its  presence.  This  geometry  is  a  typical  2-D  eddy  current  problem  and  the  plot  of 
the  induced  eddy  current  is  shown  below. 

Fig.  5  (a)  shows  the  induced  eddy  current  using  EFG  method  by  the  1st  order  polynomial 
basis  function  with  the  existence  of  the  crack  and  (b)  shows  the  image  using  EFG  method  by 
the  2nd  order  polynomial  basis  function.  Due  to  the  complexity  of  the  problem,  the  solution 
with  the  2nd  order  polynomial  basis  function  is  no  longer  convergent.  And  there  is  no  need  to 
consider  much  higher-order  basis  function. 


Figure  5.  Induced  eddy  current  plot  using  EFG  method  by  polynomial  basis  function  (a)  1st 

order  (b)  2nd  order. 
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Fig.  6  (a)  shows  the  current  plot  using  EFG  method  by  2nd  order  orthogonal  basis  and  (b) 
shows  the  difference  compared  with  Fig.  5  (a).  It  is  very  clear  that  the  2nd  order  orthogonal 
basis  function  still  converges  and  only  minor  difference  is  observed  compared  with  the  1st 
order  polynomial  basis  function. 
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Figure  6.  Induced  eddy  current  plot  using  EFG  method  by  orthogonal  basis  function  (a)  2nd 

order  (b)  comparison  with  Fig.  5  (a). 


When  using  orthogonal  basis  function,  there  is  no  need  to  perform  the  LU  decomposition 
because  of  the  avoidance  of  the  shape  function  inversion.  Fig.  7  shows  the  relative  iteration 
error  during  convergence  for  1st  order  polynomial  basis  and  2nd  order  orthogonal  basis.  We 
can  conclude  that  the  usage  of  the  orthogonal  basis  function  will  result  in  a  faster 
convergence. 


Figure  7.  the  iteration  error  versus  iteration  number  for  polynomial  basis  and  orthogonal  basis. 
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3.  Conclusion 


Conclusion:  from  the  comparison  of  the  simulation  results  in  1-D  and  2-D  electromagnetics 
problems,  the  improved  formulation  of  the  EFG  method  is  validated.  The  implementation  of 
orthogonal  basis  function  will  avoid  the  inversion  of  shape  function  matrix  and  increase  the 
stability  of  the  system.  It  provides  an  opportunity  for  the  application  of  higher-order  basis, 
which  will  increase  the  solution  accuracy  consequently.  The  above  1-D  and  2-D  examples 
have  already  shown  this  advantage  and  the  3-D  problem  will  potentially  be  simulated. 
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